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SUMMARY 


Wind shear, the variation of horizontal atmospheric winds with altitude, 
has been identified as a causative factor in several airplane accidents and may 
have been a contributing factor in others. Consequently, wind shear and its 
effect on aircraft have become the subjects of research. This study extends 
past work and concentrates on longitudinal motion. The airplane is represented 
by the three-degree-of-f reedom equations for longitudinal motion. Because 
stability and control problems occur only if the wind shear parameter 
exceeds 1 , all wind shears used in this study produced a value of greater 

than 1 . Previous investigations have shown that wind shear has little effect 
on the short-period mode for the type of airplanes used in this study. The u 
stability derivatives (derivatives with respect to perturbation velocity) were 
varied to determine the effect of changes in their magnitudes on the stability 
of the long-period mode in wind shear. It was found that increases in the 
pitching-moment derivative Mu and decreases in the vertical-force deriva- 
tive Zu made the airplane more stable in wind shear. During the study of the 
u derivatives, a wind shear tolerance factor was developed, which is a function 
of the basic stability derivatives of the airplane. If this factor is greater 
than 1, the airplane is stable in positive wind shear. 


INTRODUCTION 

Wind shear, the vertical (altitude) variation of horizontal wind, has been 
a causative factor in many airplane accidents (refs. 1 and 2) . Wind shears can 
occur at any altitude. During midcourse flight, when the airplane is at high 
altitude, an encounter with wind shear poses little danger to the airplane. 
However, during take-off and landing operations, when the airplane is close to 
the ground, an accident may occur before the airplane can recover from a wind 
shear encounter. Because low-altitude encounters represent a hazardous condi- 
tion, this investigation was restricted to the study of the problems of wind 
shear encounters on landing approach below an altitude of 130 meters. 

References 3, 4, and 5 give the results of general studies of wind shear 
for descending and climbing flight near the ground. The wind shear param- 
eter Qu, introduced in references 3 and 4 and also used in this study, is a 
measurement of the severity of the wind shear. Since a wind shear parameter 
greater than 1 produces unstable conditions, gradients for the study were 
selected to give > 1 . 

The principal thrust of this study was to examine the role the speed 
stability derivatives play in the interaction of the airplane and wind shear. 

In addition, the use of airspeed control systems to control an airplane during 
a wind shear encounter was further investigated. 

The equations of motion are derived in appendix A. A program written for 
the Hewlett Packard HP-67 programmable pocket calculator, which represents an 



improvement on program 1 of reference 6, is given in appendix B. The use of 
this calculator does not constitute an endorsement of the product by the 
National Aeronautics and Space Administration. 
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^1 f ^2/ • • 

• fSis 

coefficients of 3 x 3 matrix (eqs. (A8)) representing 
linearized equations of motion for stability; defined 
in equations (A9) 

b*| ,b2f . 

• f bg 

coefficients used in equations (A8) and defined in 
equations (A9) 
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thrust coefficient, 


PSUq2 
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3Ciji 

= , m”^-sec 

du 

c 

mean aerodynamic chord, m 

D 

d 

dt 

5 

drag , N 

Du 

1 35 

= - — , sec“* 
m 3u 

FX 

force in x-direction, N 

FZ 

force in z-direction, N 

g 

acceleration due to gravity, 9.80665 m-sec”^ 

h 

altitude, m 

Iy 

moment of inertia about y-stability axis, kg-m^ 

t + ■> 

ifDfk 

unit vectors 

i/ j/k 

indices 

/ • • • /Ky 

airspeed control system gains 

>1 

radius of gyration, m 

L 

lift, N 

M 

pitching moment, N-m 
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ly 9q 

Mu 
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ly 3u 
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V 

Va/Va 


/Vb 


1 3M 

= — — , rad“'-sec“^ 
ly 3ot 


1 8m 

= — — f rad’”' -sec*”* 
ly 90t 


1 9M 1 o 

= — f rad*"' -sec"*^ 

ly 96^ 

mass, kg 
per iod , sec 

pitching velocity, rad-sec“^ 
wing area, 

shear tolerance factor 
thrust, N 


1 9t 

= - — , sec ' 
m 9u 

time to double amplitude, sec 
time to damp to half-amplitude, sec 
forward velocity, m-sec~^ 
steady-state velocity, m-sec”^ 
perturbation velocity, m-sec“^ 

wind speed in x-direction in inertial coordinates, 

dUw 

wind shear gradient, , sec“^ 

dh 

airspeed, m-sec“^ 

airplane velocity vector and magnitude, m-sec”^ 
airplane velocity vector in body axes, m-sec~^ 
airplane acceleration vector in body axes, m-sec”^ 


m-sec ^ 
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''c 

Vr 

Vw 


W, 


w 


W, 


w 


command velocity, iti-sec"^ 
resultant velocity, m~sec“^ 
wind velocity vector, m-sec”^ 

component of in the z-direction (body axes) , m-sec"^ 

speed of updrafts or downdrafts, m-sec“^ 

dWyj 

updraft-downdraft gradient, , sec*^ 

dx 


x,y,z inertial axis system and coordinates 
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sec 
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m 9a 


Za 


1 ^^Z 

- — T“/ m-rad“ ' -sec“^ 
m 9a 


1 9^2 

Zge = ~ r m-rad l-sec~2 

m 96o 
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perturbation angle of attack, rad 
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cx^j. trim angle of attack, rad 

r total flight-path angle, Tq + /, rad 

Tq steady-state flight-path angle, rad 

Y perturbation flight-path angle, rad 

(Sg elevator deflection, rad 

throttle deflection, rad 
Cd danping ratio 

6 pitch angle, rad 

C,ri/5 moving axis system for airplane 

p air density at 0°C and 1 atm, 1.2929 kg-m"^ 

cr^ wind shear parameter for horizontal wind 

updraft-downdraft parameter 
Tg engine time constant 

03 circular frequency, rad-sec"^ 

o)j^ undanped circular frequency, rad-sec"^ 

Subscripts: 
tr tr im 

0 steady-state conditions 

Dots over a symbol indicate derivatives with respect to time. 

An arrow over a symbol indicates a vector. The symbol without the arrow 
indicates magnitude. 

AIRPLANES AND CONDITIONS OF STUDY 
Airplanes Represented in Study 

Two jet transport airplanes were represented in this study. One was a 
large four-engine, long-range jet transport referred to as airplane A. The 
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other was a small twin-engine, medium-range jet transport referred to as air- 
plane B. The aerodynamic and physical characteristics of both airplanes are 
given in table I. 

Airplanes A and B were used throughout the study. However, in order to 
give more data on the shear tolerance factor, two additional four-jet transports 
were used and are referred to as airplanes C and D. Because of the limited use 
made of airplanes C and D, their aerodynamic and physical characteristics are 
not presented. 


Wind Conditions 

Wind shear is an atmospheric phenomenon and varies with one or more atmo- 
spheric parameters. In this paper wind shear was restricted to a variation of 
the horizontal wind speed with altitude, as shown in figure 1 . If an airplane 
follows the flight path shown in figure 1 and penetrates the wind shear layer, 
then the airplane will experience a positive wind shear; that is, a head wind 
changes to a tail wind during transit of the airplane through the wind shear 
layer. Work reported in references 3 and 4 shows that unstable conditions occur 
if the wind gradient is large enough to produce a wind shear parameter greater 
than 1. For the airplanes used in this study a gradient of 0.25 sec“^ gives 
a value of a^ in the range 1.5 to 2.0. This gradient has been measured in 
thunderstorms. (See refs. 2 and 7.) 

For both stability calculations and time histories, the wind was described 
by the equation 


^w ^w,0 "*■ ^w 


(Az) 


where Az = - Zj^, with z^ being the present value of z and Zj^_i the 

previous value of z. In stability calculations taken as zero, and 

in time histories -6.10 m-sec“^ head winds are negative. For ^ 

> 0, and , this equation produces winds that provide a positive 

wind shear situation with respect to the airplane. 

In this study wind shear occurred in the altitude interval 106 ^ H ^ 56. 


Equations of Motion 

The equations of motion developed for stability calculations are derived 
in appendix A and have ^ 0 and = 0. (A program for calculating the 
characteristic equation is given in appendix B.) For the cases studied in this 
investigation, is taken as zero and equations (A4) become 
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a 

= 

26e 


Y 


«6e 

J 
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and are the equations used for the stability calculations. These equations form 
the basis of program 1 of reference 6. This program and others in reference 6 
were used to calculate all stability information. 

Time histories of spatial motions of the aircraft were generated using the 
three-degree-of-freedom nonlinear equations of longitudinal motion (appendix A, 
eqs, (A13)), The effects of wind shear, updrafts, and downdrafts were incor- 
porated in these equations, as suggested by Etkin (ref. 8). 

Gera (ref. 9) has proposed a new set of equations for the analysis of wind 
shear. These equations are restricted to the condition Fq = 0. Since the 
present investigation considered only descending flight, the new equations pro- 
posed in reference 9 were not applicable. It should be noted that for = 0 

and Tq ^ 0, the equations of this report and those given in reference 9 are 
the same , 


AIRPLANE RESPONSE TO WIND SHEAR 

Results presented in reference 4 established the effect of wind shear on 
the longitudinal stability characteristics of the airplane: Wind shear had 
little or no effect on the short-period mode, and the long-period mode changed 
from a lightly damped oscillation to two aperiodic modes, one of which became 
unstable when > 1 . A gradient of u^ = 0.25 sec“^ was used for the 
present study. The corresponding values of were 2.0 for airplane A and 

1 .7 for airplane B, Table II contains the eigenvalues and stability parameters 
of airplanes A and B for u^ = 0 (no shear) and u^ = 0.25 sec"^ . The data 
presented in table II confirm the results published in reference 4. These 
results show that for the same gradient, airplane B takes 46 percent more time 
to double amplitude than does airplane A. Material presented in figures 5 
and 6 of reference 4 indicates that the approach speed of an airplane affects 
its reaction to wind shear. Thus, the longer time required for airplane B to 
double amplitude is not unexpected, as its still-air approach speed is 
10.24 m-sec”^ less than that of airplane A, 

The wind shear data presented in table II for airplanes A and B are the 
baseline cases for determining the effects of varying the u derivatives. 

These cases are labeled "Basic airplane” in tables III and IV. 
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Effect of Changes in Magnitude of Speed Derivatives 
on Airplane Response 

The speed derivatives Dy, Z^, My, and Ty have been neglected in the 
past because these derivatives are small and contribute little to the response 
of the airplane. Because wind shear can introduce large speed changes, the 
speed derivatives were varied to determine their effect on airplane response in 
the presence of wind shear. Four values were used for each derivative: zero, 

the nominal value, 2 times the nominal value, and 4 times the nominal value. 

The speed derivatives appear in the first column of equation (1), and because 
the terms in this column as well as in the first row have little effect on the 
short-period mode, variations in the speed derivatives should affect only the 
long-period mode. Calculations showed that there was practically no effect on 
the short-period mode. Consequently, no short-period mode data are reported. 

Increasing 5y and Ty above the nominal value caused a slight increase 
in the time required for airplane A to double amplitude. (See table III.) 
Setting Dy and Ty equal to zero caused only a slight decrease in the time 
to double amplitude. Variations_in My and Zy produced more significant 
results than did variations in Dy and Ty. When the value of My was greater 
than the nominal value, the time to double amplitude increased; and at 4 times 
the nominal value of My, the aperiodic mode turned into an unstable oscillation 
with a time to double amplitude so large that control of the airplane should not 
be affected. When Zy was increased above the nominal value, the aperiodic 
mode became less stable, an opposite result compared with variations of the 
other speed derivatives. However, as Zy became smaller than the nominal 
value, the airplane became less unstable until at Zy = 0 a stable oscillation 
with an 80-sec period replaced the aperiodic mode. The process of varying the 
speed derivatives was repeated for airplane B and the results are presented in 
table IV. A comparison of the results given in tables III and IV shows that 
the trends noted for airplane A when the speed derivatives are varied are the 
same for airplane B, with the exception of the results for My. As My is 
increased, the aperiodic mode becomes a long-period stable oscillation; and at 
4 times the nominal value of My, the oscillation becomes unstable. 

The difference in response to the variations of My was traced to the 
signs of coefficients in the characteristic equation. The normalized form of 
the longitudinal characteristic equation is 


D^ + C 3 D 3 + C2d 2 + CiD + Cq = 0 


In the case of airplane A when wind shear is present, both Cg and C-| are 
negative, and increasing My causes Cg to become positive but causes C-] 
to become more negative and introduces the unstable oscillation. In the case 
of airplane B as My is increased, C-| remains positive and a stable oscil- 
lation occurs when My is twice the nominal value. When this value is doubled, 
an unstable oscillation with a time to double amplitude of about 42 sec is 
present. This difference in response to the variations in My was traced to 
the magnitude of Dy, which is larger for airplane B than for airplane A. 
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Table V shows the effect on the response of airplane A of increasing when 

= 0.0844, In this case the unstable aperiodic response becomes a damped 
oscillation at 4 times the nominal value of M^. Further increases in 
would cause the oscillation to become unstable. 

Because and have so great an influence on the airplane response 

in wind shear, a further study of these two stability derivatives was made. 

(See appendix C.) Figures 2 and 3 show the stable and unstable regions of 
and Z^ as a function of the wind shear parameter 

The magnitude of the speed stability derivatives, especially and Z^, 

has been shown to affect the response of the airplane in wind shear. However, 
it is not known whether the presence of wind shear influenced the results 
obtained or not. To obtain information on the influence of wind shear, the 
speed stability derivatives were varied in the same manner as before but with 
0u = CJw = 0. The results on the short-period mode were the same as before: 
there was practically no effect on this mode. The results for the long-period 
mode of airplane A are presented in table VI. A comparison of tables III 
and VI clearly indicates that the responses obtained with wind shear present 
are due to the interaction of the wind shear and the changes in the speed 
derivatives. Similar results were obtained for airplane B. 


Shear Tolerance Factor 

The results presented so far show that the magnitudes of and Z^ can 

have an important effect on the response of the airplane in wind shear. A 
relationship between these derivatives and and Zq^, called the wind shear 

tolerance factor Sijip, has been developed. (See appendix C.) The shear toler- 
ance factor is dependent only on stability derivatives and is 


Stf 


M^Za 


If 0 < Sijip < 1 , the airplane is unstable in wind shear and the degree of insta- 
bility becomes less as S»jip ^ ^ . If Srpp > 1 , the airplane is stable in wind 
shear . 

The emphasis has been placed on varying the u derivatives, with the most 
stress placed on My and Zy. As can be seen from the definition of Srpp, both 
Mot and Zq^ can be varied to achieve the same result as obtained from varying 
My and Zy; however, varying and Zot causes changes to occur in the short- 

period mode, which is practically unaffected by changes in Zy and My. 

The shear tolerance factor and the time to double amplitude are presented 
in the following table for four airplanes and a wind gradient of 0.25 sec“^ : 


11 



Time to double amplitude 


Airplane 

Shear tolerance factor 

A 

0.266 

C 

.481 

B 

.538 

D 

.957 


5.51 

6.59 

8.06 

28.78 


Airplanes C and D were included to extend the calculations with the shear toler- 
ance factor. Both are four-engine commercial jet airliners. Airplane C is 
medium range; p is long range. 


Spatial Motions 

Typical controls-f ixed spatial motions for a, T, h, and V for an 
encounter with a strong wind shear (a^ = 2.0) are presented in figure 4. These 
motions were generated by integrating the longitudinal equations of motion (A13) 
on a large digital computer. As can be seen, the principal characteristics of 
the motion are an increasing angle of attack and a decreasing velocity. Since 
r is related to 6 and ot through F = 0 - a, the pitch angle 9 is also 
decreasing. These conponent motions combine to cause a rapid loss of altitude, 
and for the case shown, the airplane would impact the ground before recovery 
to an equilibrium flight condition could take place. 


Control in Wind Shear 

The results presented in reference 4 show that two control systems, atti- 
tude and airspeed, are required for adequate control during wind shear encoun- 
ters. This was confirmed by results presented in reference 9. In a wind shear 
encounter (fig. 4), the angle of attack ot increases as the airspeed V 
decreases. This type of variation of a and V and an inspection of the eigen- 
vectors suggest that airspeed should be controlled. Airspeed may be controlled 
by adjusting the engine throttles, drag-only speed brakes, wing flaps, and 
finally spoilers. The foregoing systems may be used alone or in combinations. 
Reference 10 discusses and illustrates control systems using throttle position, 
flaps, and drag-only speed brakes. To evaluate the speed control as a means 
of smoothing the response of the airplane to wind shear, a control system very 
similar to the analog throttle control system discussed in reference 4 was used. 
Figure 5 is a block diagram of the airspeed control system used in this study. 

The results are presented in figure 6 for several different values of the engine 
time constants Tg, a constant that characterizes the time delay between calling 
for a given thrust level and adhering to the desired level of thrust. Three 
values of the engine time constant were used: 0.125, which would correspond 

to thrust modulators (ref. 4); 2.50, which is a good approximation for the 
engines used on airplanes A and B; and 3.75, a larger engine time constant than 
would normally be encountered. For the smallest value of the engine time con- 
stant, the response is the same as one would expect for an undisturbed airplane. 
As the engine time constant increases, the response becomes more oscillatory and 


12 



I 


less acceptable, as shown in figure 6. Although the response deteriorates, it 
does not become aperiodic up to the maximum value used for the engine time 
constant. 


CONCLUDING REMARKS 


An analytical study has been made of the longitudinal response of jet 
transport airplanes to vertical variation of the horizontal winds. The wind 
variation occurred between altitudes of 106 and 56 meters and changed a head 
wind into a tail wind (positive wind shear) . The flight condition was landing 
approach. 

When an airplane encounters wind shear, the ensuing response is a result 
of the interaction between the airplane and its environment. This interaction 
was investigated by varying the so-called u, or speed, stability derivatives. 
These derivatives were selected for variation because they affect the long- 
period mode, the mode that becomes unstable in wind shear, and have little 
effect on the short-period mode. It was found that increasing the pitching- 
moment derivative M^ or decreasing the vertical-force derivative Z^ could 
change an unstable flight situation to a stable one. In order for changes in 
My to be effective in restoring stability, it was found that the coefficient 
of the first-power term of the characteristic equation must be positive. This 
coefficient can be maintained greater than zero by increases in the drag 
derivative Dy. Alternatively, a decrease in Zy as My increases will 
achieve the same result. 

Another result of the investigation of the effects of varying the u 
stability derivatives was the development of the shear tolerance factor, which 
is a function of My, Zy, and Z^f the derivative of vertical force with 
respect to angle of attack a. If the shear tolerance factor is less than or 
equal to 1, the airplane will be unstable in positive wind shear; if the shear 
tolerance factor is greater than 1, the airplane will be stable in positive wind 
shear. The wind shear tolerance factor can be used to monitor changes in My 
and Zy and thus assure the attainment of stable conditions. 

Control of airplanes in wind is an important factor in wind shear research. 
Previous studies have reported that a flight-path control system or a pitch- 
attitude control system, as well as an airspeed control system was necessary 
for adequate control in wind shear. Results presented in this paper show that 
adequate control can be achieved through the use of an airspeed control system, 
provided the engine time constant is small enough. No other control system is 
required. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 14, 1980 
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APPENDIX A 


EQUATIONS OF MOTION 

The axes systems for the equations derived in this report are shown in 
sketches (a) and (b) . 



Sketch (a) 


The inertial axes x^y^z are shown in sketch (a); a horizontal wind 
is flowing in the positive x-direction and a vertical wind is flowing in 

the positive z-direction. Instead of being referenced to an inertial reference 
frame, the airplane is referenced to the air-mass axes. The moving axes used 
for the airplane are the system shown in sketch (b) . These are wind 

as well as stability axes for longitudinal motions, since the ?-axis is 
aligned with the single component airplane velocity. As shown in sketch (b) , 
r is measured from the local horizontal to U. As was done in references 4 
and 9, it is assumed that shear effects in the moment equation are negligible 
and that T =0 - a. 
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Sketch (b) 


The inertial acceleration expressed in the moving-axis system C/Hrt is 



+ V, 


w) 


+ 


-V > 

where = iU + j (0) + 

The wind velocities 


^ ( Va + V„) 

k(0), V„ = iu„ + j(0) + kw„ 
and w„ were taken as 


y 

and oj 


(A1 ) 


i(0) + jf + k(0) . 


“w = ^w,0 + u„(z) 


and 


«w = Ww,0 + w„(x) 


There is no loss of generality if u„^0 ''w,0 set equal to zero and 

the wind velocities taken as 


u„ = u„(z) 
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and 


“* ( ^ ) 


Thus the wind vector can be written as 


^ <3u„ ^ dw„ 

V„ = i z + k X 

dz dx 


and the acceleration due to variation in the wind speed is 


Vw = u„z + w^,x 


where 


I 



du„ 

dz 


and 


Ww 


dw„ 

dx 


and 


z = U sin r + w^ 

X = U cos r + 

-y ~y ~y- 

If the definitions of and are substituted into the expres- 

sion for the acceleration (eq. (Al ) ) , the nonlinear equations can be written as 
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U + COS r - u;u sin r cos r - w^U cos T sin T 


• • . P . p X 

- WwUy, sin ^ + g sin 1 = — 

m 


sin r - u^U sin^ V + w^U cos^ T + cos T 


- Tu - g cos r = — 
m 


' ^ 


V. (A2) 




The following assuir 5 >tions were used in the linearization of equations (A2) : 


r = To + Y 
U = Uq + U 
Y = 9 - a 


Here, 9 and a represent either the total angle of attack and pitch angle or 
perturbations from a trim condition. The latter use occurs only in linear 
equations such as (A8) . 

The wind shear parameter was also introduced. This parameter is a dimen- 
sionless acceleration whose magnitude is proportional to the wind gradient. The 

parameter = , and = . a value greater than 1 for CJ^ or 

g g 

would represent a severe wind shear, as unstable roots appear in the long-period 
mode. The linearization process given in reference 8 was followed for the 
linearization or equations (A2) . The resulting linear equations of motion are^ 


^The variables x 
transform into 5 and 
to zero. 


and 2 that appear in the definition of u^ and w^ 
C; and since xq = zq~q, both and Cq equal 
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d2 M 

— [y + a) = — 
dt2 ly 


The steady-state parts of equations (A3) are 


^X,0 

-gC?Qi sin Tq cos Fq + g sin Fq = 

m 


^Z,0 

-gdij sin^ Fq - - g cos Fg = 


\ (A4) 


0 


Mg 

ly 


J 


These last equations are used to compute the trim values of a and the thrust. 

Equations (A3) and (A4) are complete except for the definition of aero- 
dynamic forces and moments. In developing these forces and moments it should 
be remembered that they are functions of the speed of the airplane with respect 
to the air mass, so that U = Uq + u. Here Uq is the still-air approach 
speed, and u is the small speed deviation from Uq* The lift force will be 
worked out in detail. Before substitution of Uq + u for U, the force F 2 
has the form 


1 8 


PSU^ 

(Cl,o + Ci^a + Ci^a + CL^q + 
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The terms in parentheses represent the lift coefficient of the airplane and, 
except for lift coefficient at a = 0) , represent the effect of 

angular motions on the lift. The effect of speed changes, since it is con- 
tained in the dynamic pressure. When Uq + u is substituted for U and 
+ a is substituted for ot, the expression for lift becomes 

P . PS(Uq2 + 2Uqu + u2) 

Cl,o + CL(j(«tr + a) + + CLgefieJ ^ 


where is considered a second-order term and, therefore, can be neglected. 

PSUq2 

The term (Cl o *** ^i^^tr) lift required to trim the airplane, 

V ' ^ ^ / . PSUq2 

Fz,0 equations (A4) , and + Cj^q + CLgg’^ej — ^ 

lift due to angular motions from the trim state. The effect of a change in 
forward speed is given by ClPSUqu. Thus, the total lift force F 2 for equa- 
tions (A3) is 


PSUq^ . PSUq^ PSUq^ 


Fz = -[clPSUqu + Cj^a — — + Cr^d — — + Cj_q + Ct.^ 6 


'V 


PSUq^ 

Lfie^e — 


(A5) 


The longitudinal aerodynamic forces and the pitching moment can be treated 
in a similar manner. Thus, 


Fx = -|c5PSUo + C5^a 


PSUq‘ 


+ Cnr 5 


PSUq' 


D6e"e ] + Tuu 


(A6) 


where 


“ l^D,o 


p2 

CL,tr 


TTA 


+ 


and Ty is the change in engine thrust due to a change of forward speed. The 
aerodynamic pitching moment is given by 
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PScUq^ 


” 1 Cj^PScUgU + C™ a 


. PScUo^ 


PScUq^ 


2 * T 


where 


j^m,o + ■*■ 




and Cni,o is- the pitching -moment coefficient for a = 0. Equations (A5) , (A6) , 
and (A8) are substituted into equations (A3) to obtain the linearized equations 
of motion that are used for stability studies. Note that in this substitution 

u = w = — / and a = ^ ; the equations then take the form 
dt dt Uq dt 


+ aiD + a2 


C] r bi b2 


agD^ + a 7 D + as agD + a-| q ? = t>3 t>4 (AS) 

ai 2D^ + ai 3D^ + ai 4D + 3150 Y bs bg 


( g^T CdPSUq Ct^PSUq 

a-] = - I sin 2 Fq - + 




CBaPSUo^ 


(A9a) 


(A9b) 


(A9c) 


34 = g(cos Fq - Prp cos 2 Fq) 


(A9d) 


9 / . 2 r N 

as . - -(Ot sin2 r„ - 0„) - — 


(A9e) 
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> 

PSUq2 



^2mUo 

(A9f) 

PSUq2 

2mUo 


(A9g) 


(A9h) 


39 = - Uq - Cr 


*3 2m 


aiO = g(sin Fq - a<p sin 2Fo) 




4 “ *-n 


»15 = ■ Cmq 


4. r* • 1 

PScUq^ 

[ "tetj 

' 2IyUo 

PScUq2 


2IyUo 


PScUq^ 


2IyUo 



(A9n) 


(A9o) 


■ ■ “"“Se liT 


(A9p) 
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PSUq" 




PSUo“ 
® 2m 


b4 = 0 


^5 = Cj^e 


PScUq2 

2ly 


*^6 - 


PScUq' 

2Xv 


(A9q) 


(A9r) 


(A9s) 


(A9t) 


(A9u) 


If these equations reduce to the usual form of the linear longi- 

tudinal equations of motion. 

The characteristic equation of the system described by equations (A8) is 
obtained by expanding the determinant of the 3x3 matrix that occurs on the 
left side of this equation. The general form of the characteristic equation 
is 


CgD® + C 5 D 5 + C 4 D 4 + C 3 D 3 + C 2 D 2 + C^D + Cq = 0 (AlO) 


The constant term Cq is equal to 0^ so that equation (AlO) is a sixth-degree 
equation with a zero root. Equation (AlO) was written 


CgD^ + C5D4 + C4D3 + C3D2 + Cl = 0 


(All) 


for programming purposes. The expressions for the coefficients of the charac- 
teristic equation are 


^6 = - ai 2 ag (A12a) 

C 5 = a] (ag - a-|2a9) + (^7 ^15^6 “ ^13^9 ~ 310^12) (A12b) 
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C 4 = a-| (37 + a-jsag - a-j 3 a 9 - aioai2) + (^8 ^IS^? “ si4®9 ” ®133l0) 

+ a2(ag - a-j2a9) - a335 + a^aiz^i2 (A12c) 

C3 = a-j (ag + a-|5a7 - a-|4a9 - a-j3a-io) + ” ®14^10) 

+ 32 (a 7 + a-jgag - 3-1339 - aioa-|2) - 33(31535 - 31139) 


+ 34(31335 - aiiag) (AT 2 d) 

C2 = 31(31538 - ai 43 io) + 32(30 + 31537 - 31439 - 313310) 

+ 33311310 + 34(35314 - 31137) (A 12 e) 

Cl = 32(31530 - 3 i 43 io) - a 4 aii 38 (A 12 f) 

Co = 0 (A 12 g) 


If Oy and/or a„ = 0, the Ci coefficient becomes 0, and the charac- 
teristic equation reduces to a quartic. In this case the C 2 coefficient 
corresponds to the Co coefficient obtained from equations ( 1 ). 

An HP-67 program that calculates the coefficients of the 3 x 3 matrix 
(eqs. (A 8 )) and the coefficients of the characteristic equation (eqs. (AlO) 
and (A11 ) ) is given in appendix B. 

The equations of motion presented in this appendix model the condition in 
which an airplane is subjected to a combination of vertical wind shear, updrafts, 

g^OuOw^ g^o^OwC 

and downdrafts and contain wind interaction terms and not 

present in the previous formulation of this problem. If either or is 

set to zero, these equations reduce to equations given in reference 4 when the 
same parameter is set equal to zero. If the wind interaction terms are neglected 
and and are not equal to zero, the resulting equations will predict a 

larger instability (a shorter time to double amplitude) than equations (A 8 ) . 

Equations (A 8 ) are approximate rather than exact equations. The inexact- 
ness occurs because of the assumption that F = 6 - a. The angle T, defined 
in sketch (b) , is not the traditional F because it is not inertial. The 
traditional F, which will be called Fj, can be obtained by applying the cor- 
rection Fj = F - Sg, where Sq is given by 

( -Uy^ sin F - Wy, cos F 

— 

U - u^ cos F + w^ sin 
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Since £g can be considered a small angle, it can be set equal to the 
argument of the inverse tangent. This correction and its first and second 
derivatives would be applied to equations (A3). However, equations (A3) were 
used with equal to zero; thus, the correction would be 


u^, sin r 
U - U;, cos r 


which approaches zero as F approaches zero. The corrections for the first 

r 

and second derivatives are proportional to — and are smaller corrections 

u2 

than Gg* Because the corrections and F are small, no corrections were 
applied when equations (A8) were used to calculate the stability parameters. 

If updrafts and downdrafts had been included in the analysis, corrections to 
r would have been necessary. Equations (A8) were not used to compute time 
histories. 

The six-degree-of-f reedom equations of motion for principal body axes are 
m^Vg + 00 Vg y — F 

H + 00 X H = M 

where H = IxxP ^ ^ total force and moment vectors; 

p, q, r are the angular velocities; and are the products of inertia. 

These equations may be separated into lateral and longitudinal equations of 
motion; the latter are 


^x 

u + qw = - g sin 0 + — 

m 


w - qu = g cos 0 + — 

m 


0 = M/Iy 


) (A13) 




where u, w, q, and 0 are total values, not perturbations and Fy^, F^r 
and M are the total aerodynamic forces and moments acting on the airplane. 
Equations (A13) were used to compute the time histories presented in this 
report. Wind shear was accounted for by making the wind component of the 
airspeed vector altitude dependent. 
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PROGRAM FOR CALCULATING THE CHARACTERISTIC EQUATION 
OP LONGITUDINAL MOTION 

This program calculates the coefficients of the 3x3 matrix in equa- 
tions (A8) using the expressions given in equations (A9) . After the calcu- 
lation of these coefficients has been completed, equations (A12) are used to 
calculate the coefficients of the characteristic equation (AlO) under certain 
circumstances. The characteristic equation may be fifth order, so that a 
routine to extract the real root of a quintic is included in the program. 
Programs 4 and 5 of reference 6 can be used to complete the calculation of the 
roots of the characteristic equation and to determine the stability parameters. 

The two-card program that follows is an improved version of program 1 of 
reference 6. This program may be checked by using the checks case for pro- 
gram 1 of reference 6. 
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Card 1 


Step Key Entry Comments 

001 LBLA 

RCL8 
STO+7 
STOx9 
RCLO 

CHS 

STO-r(i) 

RCL7 

01 0 RCL6 

SIN 
x^ 

X 

RCL8 

STOO 

RCL6 

COS 

ST08 

020 + 

RCL4 

X 

RCLB 

RCL2 

X 

RCL5 

X 

RCLl 

030 STOE 

RCL5 

X 

2 

• 

STOx3 

STOB 

T Cl 

p>s 

ST05 

040 RCLD 

RCLl 

f 

RCL6 

X 

STO+9 


Step Key Entry Comments 

RCL5 

x2 

RCLC 

050 T 

IT 

STO+4 

RCL2 

STO+3 

RCL8 

STO+7 

RCLE 

STOx4 

060 STOX5 

RCLA 

X 

STOx9 
RCLB 
STOx (i) 

STOxO 

STOxl 

STOX2 

STOX3 

070 RCLA 

STOx6 
STOx7 
STOx8 
P»-S 
RCL8 
RCL6 
2 
X 

COS 

080 RCL7 

STOx8 

X 

STOA 

RCL6 

SIN 

STOx8 

RCL8 

2 

090 X 

RCL4 
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Step Key Entry Comments 
STOx (i) 

X 

STOB 

RCL4 

RCL5 

STOX8 

TOO STOxO 

STOX9 
RCL9 

X 

STOC 

RCL3 

RCL8 

+ 

RCLO 

RCL5 

no p>- S 

CHS 
STO+2 
CHS 
STCH-0 
STOfI 
STOf3 
STOf6 
STOt 7 
1/x 

1 20 p+S 

ST06 
I«-S 
R'l' 

STO-5 

Ri 

STO-4 

1 

0 

STOI 

1 30 RCL3 

P»-S 
RCL6 
P>-S 
RCL2 

X 

STOD 

P>S 

ST05 C6 

1 40 p>S 

RCL4 

X 


Step Key Entry Comments 

RCLl 

RCL8 

RCL3 

X 

+ 

RCL2 

RCL7 

1 50 X 

RCLB 

P>S 

RCL6 

P>S 

X 

STOE 

+ 

160 STO(i) C5 

ISZ 
RCLD 
RCLC 

X 

RCLE 

RCL4 

X 

+ 

RCLO 

1 70 RCL5 

X 

RCLA 

RCL5 

P>S 

RCL6 

P>S 

X 

X 

1 80 + 

RCLC 

RCL8 

RCLl 

X 

+ 

RCL6 

RCL2 

X 

1 90 RCL7 

RCLB 

X 
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Step 


Key Entry 


STOD 

+ 

STO(i) 

RTN 


Comments 


C4 
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Storage Map for Card 1 


Address 

Register 

Initial storage 

Output storage 

0 

Rq 

ky 

^5 

1 


m 

b4 

2 

R2 

P 


3 

R3 

Cni 

■'•u 


4 

R4 

g 


5 

Rs 

Uo 

be 

6 

Re 

To 


7 

R7 



8 

Rs 

<^w 


9 

Rg 



10 

So 

^Dot 

33 

n 

Si 


a? 

12 

S2 

Cl0 

ag 

13 

S3 


ae 

14 

S 4 

0 

pi 

0 

ai 

15 

S5 

0 

as 

16 

Se 

^ma 

ai4 

17 

S7 

^i"d 

ai3 

18 

Ss 

‘-me 

a-15 

19 

Sg 

*-m,o 

an 

20 

Ra 

c 

a4 

21 

Rb 

S 

aio 

22 

Rc 

A 

^2' ag 

23 

Rd 

0 

u 


24 

Re 

0 


25 

I 

20 

ai2 
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Card 2 


Step Key Entry Coniments 

001 LBLA 

ISZ 
RCLC 
RCLE 

X 

RCLD 

RCL4 

X 

+ 

01 0 RCL 8 

RCLC 

X 

RCL 6 

RCLB 

X 

STOE 

+ 

RCL 8 

020 RCL5 

X 

RCL9 

RCL2 

X 

RCLO 

X 

RCL7 

030 RCL5 

X 

RCL9 

RCL3 

X 

RCLA 

X 

+ 

STO(i) C 3 

040 ISZ 

RCLD 
RCLC 

X 

RCLE 

RCL4 


Step Key Entry Comments 
+ 

RCLO 

RCL9 

050 RCLB 

X 

X 

+ 

RCL5 

RCL 6 

X 

RCL9 

RCLl 

X 

060 

RCLA 

X 

+ 

STO(i) C 2 

RCLC 

RCLE 

X 

RCLA 

RCL9 

070 X 

RCLC 

X 

P>S 

ST04 Cl 

RCL5 

STOtO 

STOtI 

STOt2 

080 STOt3 

STOt 4 
RCL4 
x =0 
GOTOC 
GOTOD 
LBLC 
4 

RTN 

LBLD 

090 5 

RTN 

LBLB 


X 
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Step Key Entry Comments 

RCLO 

STOA 

RCLl 

STOB 

RCL2 

STCX: 

RCL3 

1 00 STOD 

RCL4 
STOE 
0 

ST07 

FIX 

RCLA 

5 
• 

no 0 

8 

ST05 

ST06 

GSBb 

ST08 

• 

1 

6 

1 20 STO+6 

LBLO 
1 

STO+7 

GSBb 

ST09 

RND 

Pause 

x=0 

G0T01 

1 30 RCL6 

RCL5 
RCL6 
ST05 
x^y 

RCL8 

RCL9 

ST08 

x>y 

1 40 


Step Key Entry Comments 

X 

ST06 

GOTOO 

LBLl 

RCL7 

• 

1 50 6 

9 
3 

CHS 

RCL6 

• 

RCL6 

5 

RTN 

LBLb 

1 60 2 

0 

STOI 

1 

GSBc 

STOO 

GSBc 

STOI 

GSBc 

ST02 

1 70 GSBc 

ST03 

GSBc 

ST04 

RTN 

LBLc 

RCL6 

X 

RCL(i) 

+ 

1 80 ISZ 

RTN 


RCL8 
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Storage Map for Card 2 


Address 

Register 

Initial storage 

Output 

Output 


LBLA^ 

LBLB 

0 

Ro 

The initial 

C 5 

d3 

1 

Ri 

storage for 

C4 

<32 

2 

R2 

card 2 LBLA 

C 3 


3 

R3 

is the 

C2 

do 

4 

R 4 

output of 



5 

R5 

card 1 . 



6 

Re 




7 

R? 




8 

Rs 




9 

Rg 




10 

Sq 


^3 


11 

Si 


a? 


12 

S2 


ag 


13 

S3 


ae 


14 

S 4 




15 

S 5 


^5 


16 

Se 


a-, 4 


17 

S? 


ai3 


18 

S8 


^15 


19 

S9 


an 


20 

Ra 


34 


21 

Rb 


aiO 


22 

Rc 


^ 2 ' as 


23 

Rd 




24 

Rg 


ai2 


25 

I 





^At end of label A, the numeral 4 or 5 is displayed in x; 5 indicates 
a quintic and 4 a quartic. If 5 is displayed, push B; if 4, insert quartic 
program and proceed. 

Note: In addition to the stored output, the stack contains the 

following information: 

Stack register 

T number of iterations to obtain root 
z time to halve or double amplitude 
y root of quintic equation 

X 5 indicates that a real root of a fifth-order equation 
has been found 
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EXPRESSIONS FOR Cq = 0 

The nonnormalized characteristic equation of longitudinal motion is 

C 4 D 4 + C 3 D 3 + C2d 2 + C-iD + Co = 0 (Cl ) 


When wind shear is introduced in the airplane stability and control problem, 
the coefficients affected most are C-; and Cq. The coefficient Cq is the 
dominant one. For negative wind shear, the airplane remains stable; however, 
for positive wind shear, the airplane becomes unstable when the wind shear 
parameter exceeds 1. (See ref. 4.) The expression for Cq (when = 0) 

is 




(C2) 


If Cq is set equal to zero, equation (C2) may be solved for and M^. The 

expression for 2^;^ that satisfies Cq = 0 is 


34^7311 " ^10(^3^11 ■ 3l3l4> 

Zu ______ + 

34811 

and a similar expression can be obtained for 
make figure 2. 

If equation (C2) is solved for the ratio 
Sion results: 



sin^ Tq 


(C3) 


Equation (C3) was used to 
ail ^"kx 

= , the following expres- 

®14 


34 

3i 35 

an aio 

= (C4) 

314 ^4 

33 37 

aiO 
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Equation (C4) may be written in the form 


= 1 (C5) 


If the equality of equation (C5) is fulfilled, Cq = 0. If the left-hand side 
of equation (C5) is less than zero, unstable conditions exist; if greater than 
zero, the airplane is stable. The matrix elements a-| and a 3 may be a 
negligible contribution to the term in parentheses of equation (C5) . If a-| 

and a 3 are dropped, the ratio in the parenthetical term cancels, and 

^10 

equation (C5) may be written as 


a-|i a-j 

« 1 (C6) 

®14 ®5 


If ai and a 3 are retained, the value in parentheses for airplane A is 
0.00466; if dropped, the value is 0.00463. This small difference justifies 
dropping a-) and a 3 . Except for a 5 , all the a-terms in equation (C6) are 
functions solely of the stability derivatives. Wind shear is thus introduced 

g - _ 

through the inclusion of the term - — o<i> sin^ To - <^w in Zu. For reason- 

Uo 

able values of Fq (that is, -0.17453 ^ Fq ^ 0.17453) and = 0, the 

expression - — sin^ Fq - will make a negligible contribution to a 5 

because of sin^ Fq in the expression. For airplane A with = 0, 

Fq = -0.05236 rad, and Oy = 2,0, this term has a value of -0.0007; and since 
Zu is 0.2553665, the wind shear term may be safely dropped when cf^ = 0. If 
Oyf is nonzero, the wind shear term for is <7 ^/Uq and, unless is 

an a7 

very small, may not be neglectable with respect to Zy. The term — is 

ai4 as 

called the wind shear tolerance factor S>pp and may be written as 
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*> 1 

%as 


or 


Mu^ 

« 1 



The interpretation of the wind shear tolerance factor is discussed in the text. 
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TABLE I.- PHYSICAL AND AERODYNAMIC CHARACTERISTICS 


FOR AIRPLANES USED IN THIS STUDY 

A constant air density of 1 ,2929 was 

L used for all calculations 


Characteristic 

kyr m 

c, m 

S, 

A 


m , kg 

g r iri'^sec^ • , , • , 

^u 



Uq, m-sec"^ . . . . 

r, rad 


cBct' 

rad"'’ 

CLa' 

rad"’’ 

•^Le 

. . . 


rad"'’ 

Cd,o 

. . . 


rad“^ 

•"me 

. . . 


Airplane A 
10.463784 
7.0104 
267.1 
7.03 
90909.1 
9.80665 
0 to 2.0 
0 

77.12 

-0.05236 

0.529 

4.87 

0.283 

0.0889 

0.038 

-0.241 

-0.707 


Airplane B 
5.38 
3.41 
91 .04 
8.83 
40916-87 
9.80665 
0 to 1.7049 
0 

66.88 

-0.05236 

0-9168 

6.59 

0.2208 

0.20893 

0.13778 

-0.0656 

-0.6884 


Cma' cad“‘> 

^m,o 

Ct^/ m"’^-sec . . . . 

^D6e' . . . . 


-1 .115 
-0.015 
-0.00025 
0.23302 
0.01387 
-1.017 


-1.369 

-0.14 

-0.0055409 

0.46413 


-1.656 



TABLE II.- RESPONSE OP 

[ro ^ 


Mode Root ^ 1 / 2 ' 

Airplane A, 

Short-period -0.6990 ± 0.8061 i 0.991 

Long-period -0.00648 ± 0.1291 i 106.97 

Short-period -0.6933 ± 0.7939i .995 

Long-period -0.1373, 0.1258 

Airplane B, 

Short-period -0.6082 ± 0.9882i 1.14 

Long-period -0.0298 ± 0.1202i 23.33 

Short-period -0.6074 ± 0.972U 1.14 

Long-period -0.1352, 0.0863 








TABLE III.- EFFECT OF VARYING THE U-STABILITY DERIVATIVES ON THE 


LONG-PERIOD MODE OF AIRPLANE A 
[Uq = 77.12 m-sec~^ ; = 0.25 sec”^ ; Fq = -0.0524] 


Parameter varied 


Root 

t■^/ 2 , sec 

tjj, sec 

P, sec 

tOjj, rad-sec“^ 


Remarks 


Du = 0 
Du = 0.0212 
Du = 0.0424 
Du = 0.0848 


Effect of varying the magnitude of D^ 


•0.1253, 0.1354 
-0.1373, 0.1258 
•0.1502, 0.1171 
-0.1784, 0.1021 



Basic airplane 


Effect of varying the magnitude of T^. 


Tu = 
Tu = 
Tu = 
Tu = 


-0.00281 

-0.00562 

-0.01124 


-0.1356, 

0.1271 


5.45 




-0.1373, 

0.1258 

5.51 



-0.1389, 

0.1246 

5.56 



-0.1423, 

0.1223 
1 

5.67 

1 



Basic airplane 


Effect of varying the magnitude of Mu 


Mu = 0 

Mu = 0.000997 
Mu = 0.00199 
Mu = 0.00399 


-0.1666, 0.1425 
-0.1373, 0.1258 
-0.1040, 0.1050 
0.0125 ± 0.0356i 


4.86 

5.51 

6.60 

55.35 


176.67 


0.0377 


-0.31 


Basic airplane 


Effect of varying the magnitude of Z, 


Zu = 
Zu = 
Zu = 
Zn = 


0 

0.255 

0.502 

1.021 


-0.00639 ± 0.0782i 
•0.1373, 0.1258 
-0.2103, 0.1970 
-0.3186, 0.2959 


108.42 



80.34 


0.0785 


0.0815 


Basic airplane 


U) 



IQ iQ iQlQ 



TABLE IV.- EFFECT OF VARYING THE U-STABILITY DERIVATIVES ON 
LONG-PERIOD MODE OF AIRPLANE B 
[Uq = 66.88 m-sec“^ ; u^ = 0.25 sec“^ ; Tq = -0.0524 rad] 


Parameter varied 

Root ^1/2^ sec rad-sec"^ Cd 

1 Remarks 


Effect of varying the magnitude of D 


0 -0.1090, 0.1023 
0.0427 -0.1352, 0.0863 
0.0955 -0.1662, 0.0740 
0.1709 -0.2365, 0.0573 


6.79 

8.06 

9.36 

12.16 


Basic airplane 


0 

0.0039 

0.0078 

0.1425 


Effect of varying the magnitude of 




0 

0.0356 

0.0713 

0.1426 

-0.1122, 0.0993 
-0.1352, 0.0863 
-0.1608, 0.0758 
:-0.2182, 0.0607 


7.00 

8.06 

9.12 

11 ifi 





1 1 • D 
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TABLE V.- LONGITUDINAL RESPONSE, AIRPLANE A 
[By = 0.0844] 
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TABLE VI.- EFFECT OF VARYING THE U-STABILITY DERIVATIVES ON 


LONG-PERIOD MODE OF AIRPLANE A 


[Uq = 77.12 m-sec~l ; = 0.0 sec~^ ; Tq = -0.0524 rad] 


Parameter varied 
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(i)u» rad-sec ^ 


Remarks 


Effect of varying the magnitude 

of Du 



Du = 0 

-0.0039 ± 0.1288i 
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Figure 1 Wind shear. 
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(a) Airplane A. 

Figure 2.- Curves for Cq = 0 for airplanes A and B as a function 

of 0,, and M,,. 
















Figure 5.- Block diagram of airspeed control system. 
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